Enhancement of kinetic energy fluctuations due to expansion 
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Global equilibrium fragmentation inside a freeze out constraining volume is a working hypothesis 
widely used in nuclear fragmentation statistical models. In the framework of classical Lennard Jones 
molecular dynamics, we study how the relaxation of the fixed volume constraint affects the posterior 
evolution of microscopic correlations, and how a non-confined fragmentation scenario is established. 
A study of the dynamical evolution of the relative kinetic energy fluctuations was also performed. 
We found that asymptotic measurements of such observable can be related to the number of decaying 
channels available to the system at fragmentation time. 
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I. INTRODUCTION 

The multifragmentation phenomenon has attracted the 
attention of the intermediate/high energy nuclear com- 
munity since more than a decade. Starting in the 1980's, 
a lot of attention has been paid to the identification and 
characterization of the occurrence of a nuclear liquid-gas 
phase transition. On one hand, indications of critical be- 
havior were obtained assuming a Fisher-like scaling re- 
lation for fragment mass distributions (il yi On 
the other, it has been claimed that the observed critical 
behavior in static observables is compatible with a first 
order phase transition 0, IE HIE 13 ■ 

As in nuclear fragmentation experiments one only 
has access to asymptotic observables, the fragmentation 
stage of the dynamic evolution must be reconstructed 
using simplified models. To accomplish this task, the 
most widely used statistical models of the MMMS [T^ or 
SMM Tl'l type, replace the intractable problem of quan- 
tal N-body correlations with a non interacting gas of frag- 
ments confined with sharp boundary conditions |lOl . 
A distinctive feature of these equilibrium fragmentation 
models is that a vapor-branch can be recognized in the 
system caloric curve (CC). Moreover, in the statistical 
sampling of configurations, the observed asymptotic flux 
is exclusively associated with the mean Coulomb energy 
of partitions inside the constraining volume. Collective 
motions are completely disregarded before and during the 
fragment formation stage. 

On the contrary, when fragmentation is considered in 
a non-constrained free-to-expand scenario qualitatively 
different features are found. In this case a local equilib- 
rium picture replaces the global one, and a characteris- 
tic flattening of the caloric curve at high energy can be 
found 12 • This happens as a direct consequence of the 
presence of a collective degree of freedom: the expansive 



motion behaves like a heat sink, and precludes the disor- 
dered kinetic energy associated with the system tempera- 
ture to increase without bounds ^El ■ These last remarks 
rise some reasonable concerns regarding the sharp con- 
straining volume working hypothesis that are worth to 
be analyzed. 

In recent contributions the study of the behavior of 
the relative kinetic energy fluctuations, Ak = Na\/ < 
K >^, has received considerable attention. It has been 
shown 0,01 that partial energy fluctuations can be used 
to measure the EOS, and to extract information of the 
possible phase transition occurring in finite systems like 
multifragmenting nuclei. Negative specific heats are ex- 
pected to signal first order phase transitions in these sys- 
tems IE ^3 ■ Taking into consideration that kinetic en- 
ergy fluctuations and the system's specific heat are re- 
lated by 0: 
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(where /3 = 1/kT) negative values of the specific heat 
are expected to appear whenever gets larger than the 
corresponding canonical expectation value: 

Using a classical molecular dynamics model it has been 
already shown in ^3 that large relative kinetic fluctua- 
tions do take place in confined classical isolated systems, 
provided that the constraining volumes are large enough 
to accommodate non-overlapping configurational frag- 
ments. In this paper we want to pursue with the system's 
characterization by studying its temporal evolution once 
we remove the constraining walls within which it was 
equilibrated. Our aim is to understand how much mem- 
ory of the initial configuration is kept in the asymptotic 
fragmentation patterns which are accessible experimen- 
tally. To that end we will focus on the dynamical evolu- 
tion of fragment mass distribution functions (MDF), and 
the relative kinetic energy fluctuation, Ak, that turns 
out to be a useful tool to probe the unconstrained frag- 
mentation scenario. 
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II. THE MODEL 

The system under study is composed by excited drops 
made up of 147 particles interacting via a 6-12 Lennard 
Jones potential with a cut-off radius = 3(t. Ener- 
gies are measured in units of the potential well (e), a 
characterizes the radius of a particle and m is its mass. 
We adopt adimensional units for energy,length and time 
such that e = o" = 1, <o = \/ cr^m./48e. Our particles 
are uncharged and the inclusion of Coulomb is certainly 
a necessary step to insure that our results are pertinent 
to the nuclear fragmentation problem However, the 
fact that the the nuclear force and the van der Waals in- 
teraction both share the same general features (i.e. short 
range repulsion plus a longer range attraction) supports 
the use of this simplified classical model to obtain quali- 
tatively meaningful results in the nuclear case, and serves 
also to explore phenomena of greater generality. The ini- 
tial condition is chosen as a spherical confining 'wall', 
introducing an external potential V^aii ^ {r ~ r^aii)^^'^ 
with a cut off distance rcut = Ic. The set of classi- 
cal equations of motion were integrated using the well 
known velocity Verlet algorithm, which preserves volume 
in phase space 19|. 

Using this system, we have performed the following 
molecular dynamics (MD) experiments. First, we let 
the system equilibrate inside the constraining volume 
at a chosen density. We call this initial condition 
the constrained state. Afterwards, we remove the 
constraining walls, and let the system evolve for a very 
long time, such that the fragment composition reaches 
stability. This will be referred to as the asymptotic stage. 

One of the key observables to study phase transforma- 
tions from a morphological point of view is the fragment 
mass distribution. The simplest and more intuitive clus- 
ter definition is based on correlations in configuration 
space: a particle i belongs to a cluster C if there is an- 
other particle j that belongs to C and |ri — rj| < rd, 
where rci is a parameter called the clusterization radius. 
In this work we took rd = rcut = 3cr. The algorithm 
that recognizes these clusters is known as the "Mini- 
mum Spanning Tree" (MST). The main drawback of 
this method is that only correlations in r-space are used. 
MST clusters give no meaningful information about the 
fragmentation dynamics when the system is dense. 

A more robust algorithm is based on the analysis of the 
so called "Most Bound Partition" (MBP) of the system, 
introduced in Ref. [23|. The MBP is the set of clusters 
{Ci} for which the sum of the fragment internal energies 
attains its minimum value: 

{a} = argmin{c.} [ii;{c,} 

JGC, 3,keCij<k 

where the first sum in (jSJ is over the clusters of the par- 
tition, Kj"^ is the kinetic energy of particle j measured 



in the center of mass frame of the cluster which contains 
particle j, and Vij stands for the inter-particle poten- 
tial. It can be shown that clusters belonging to the MBP 
are related to the most-bound density fluctuation in r-p 
space 20]. The algorithm that finds the MBP is known 
as the "Early Cluster Recognition Algorithm" (ECRA). 

In a previous contribution we studied the EOS of a 
finite, constrained, and classical system using cluster dis- 
tributions properties in phase space and in configura- 
tional space [iTj. Figure ^ summarizes our findings. In 
the upper panel we show the density dependence of the 
system energy Ecru at which transition signals (power- 
law mass spectrum, maxima in the second moment of 
the mass distribution, and normalized mean variance) are 
detected. Given a confined system with a density value 
p, for E < Ecritip) a big ECRA cluster can be found, 
whereas for E > Ecritip) a regime with high multiplic- 
ity of light ECRA clusters dominates. The lower panel 
shows the same transition line in the most usual (T, p) 
plane (see [13 for details). 

According to the calculated critical exponents r, and 7, 
three density regions could be identified. One (very low 
density limit, region labeled A in Figure^ in which the 
boundary conditions are almost ineffective, and a back- 
bending of the caloric curve can consequently be observed 
even in the isochore ensemble fJll . Another one (region C 
in Figure^ , corresponding to the high density regime, in 
which fragments in phase space display critical behavior 
of 3D-Ising universality class type. And an intermediate 
density region (region B in Figure ^ , in which power- 
laws are displayed but can not be associated to the above 
mentioned universality class (see Ref. 0|). The thermo- 
dynamic critical point of the liquid-gas phase transition 
(pc ~ . 35(7^^, Tc « 1.12eo) corresponds to the meeting 
point of these two last regions. 

Having these results in mind, we chose the follow- 
ing density values, pi = 0.026cr~^, p2 = O.lOcr"^, 
P3 — 0.6(7"^ (i.e. one from each relevant region of the 
EOS) in order to study the effect of the dynamic evo- 
lution on the observables that characterize the initially 
equilibrated confined systems. 



III. TIME EVOLUTION OF PARTICLE 
CORRELATIONS 

In order to explore the effects of constrain removal on 
the fragment spectra of the systems under study, we com- 
pare the confined ECRA-fragment mass distribution of 
the constrained system, against the corresponding MST 
asymptotic ones. 

In Figure 121 we show MDF calculations performed over 
systems initially constrained at the density values above 
mentioned (from left to right: pi,p2, and ps) with dif- 
ferent total energies (increasing E from bottom to top). 
For a given set of values {E, p) the confined state ECRA- 
fragment mass spectra and the corresponding asymptotic 
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FIG. 1: Upper panel: energy Ecru at which transition 
signals are detected in a confined, finite size fragment- 
ing system, as a function of the system density. Lower 
panel: same transition line in the temperature-density 
plane (see text for details). The arrows indicate the 
density values considered for the present calculations. 
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FIG. 2: Confined ECRA vs asymptotic MST 
mass spectra calculated for 3 different initial den- 
sities: p = 0.026cr-3,0.1cT-3, and 0.6a-^ shown 
in first, second, and third columns respectively. 



MST-fragment mass spectra can be compared. It should 
be kept in mind that for p2 and pa initial states present 
just one big fragment according to the MST algorithm. 

For the case of an initially extremely diluted system 
(first column in Figure HJ confined states are already 
fragmented in configurational space. The corresponding 
asymptotic fragment structure remains essentially unal- 
tered, aside of the evaporation of light clusters from the 
heaviest fragment. 

On the other hand, it can be seen that for initial higher 
densities {p2, and p^) the q— p correlations of the con- 
fined state, probed by ECRA partitions, are modified by 
the dynamics that follows the walls removal. For these 
densities, the size of the biggest clusters detected in both, 
the confined and asymptotic stages, seem to agree to a 
high degree. However noticeable differences are found in 
relative abundances, specially for intermediate size frag- 
ments (IMF). 

The development of an expansive flux modifies the 
correlation pattern observed in phase space for initially 
dense confined configurations acting as a memory 
loosing mechanism of their initial fragment structure. 
This result implies that even a subcritical density 
relatively diluted as the p = O.lcr^^ case cannot be 
taken as a freeze out configuration, since freeze out by 
definition implies a persistency of the information up to 
asymptotic times. It is then clear that between confined 
states and corresponding asymptotic configurations, the 
intermediate stage, characterized by the development of 
expansive collective motion and the consequent expan- 
sion of the system, plays a crucial role in the process of 
fragment formation. The relative similarity between the 
asymptotic distributions originated from the p = O.lcr"'^ 
and p = 0.6(T~^ initial conditions, suggests that the 



expansion dynamics leads to a fragmentation pattern 
essentially determined by the total energy. On the 
other hand the distribution originated from the initially 
most diluted system is sensibly different from the denser 
cases. Indeed the size of the heaviest fragment is smaller 
and the IMF distribution is steeper in this case. This 
means that a fragmentation pattern originated from an 
expansion process is not equivalent to an equilibrated 
freeze out with sharp boundary conditions. It has been 
proposed|22j| that such a fragmentation pattern may still 
be described in statistical terms through information 
theory with the constraints of a fiuctuating volume and 
a radial fiow. The possible pertinence of such a scenario 
is left to future investigations. 



IV. KINETIC ENERGY FLUCTUATIONS 

As was stated in the introduction, fragmenting finite 
systems at microcanonical equilibrium display negative 
specific heats whenever relative kinetic energy fluctua- 
tions surpass the canonical expected value. This is the- 
oretically expected in the liquid-gas phase transition of 
small systems if the volume is not constrained by sharp 
boundary conditions [2^. Consistently it has been sug- 
gested that the magnitude Ak can be used to identify 
the occurrence of first order phase transitions in such 
systems 0] • If the equiprobability of the different config- 
uration microstates is assumed, as it is by construction 
the case in our initial conditions, K = K{T), and 
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FIG. 3: Upper panels (3. a, 3.b, 3.c): MST fragment spectra calculated at t = 20,30,40 and, 
300^0 for a p = 0.6(t~^ system at three different energies. Middle panels (3.d, 3.e, 3.f): (red 
circles) and MST-biggest cluster mass (dashed lines) as a function of time. Lower panels (3.g, 3.h, 
3.i): Total kinetic energy (solid line) and Kinetic energy fluctuations (red circles) as a function of time. 



In this way, the system specific heat can be calculated 
through Eq.Q] knowing both, < K >, and crj^. 

Since the dynamics of the expansion tends to drive the 
system into a rather rarefied state, and allows event by 
event volume fluctuations which are suppressed by the 
sharp boundary, it is interesting to see whether such ab- 
normal fluctuations can be reached in a fully dynamical 
simulation. In this case the average kinetic energy can- 
not be associated to a temperature any more, and the 
relationship expressed in Eq.^can no longer be invoked 
to estimate the system specific heat. However meaning- 
ful information can still be obtained from Ak when we 
resort to dynamical quantities only. In Figure |31 we show 
the time evolution of MST mass spectra (panels|31a,|2|b, 
and|3|c), Ak (panels|nid,|3|e, andOlf), K, and crx (panels 
ElgjElli, and|31i) calculated for a system initially confined 
at a density p = — O.Ga^^, for three different energies. 
We also show in panels |21d, |3e, andElf, the correspond- 
ing time evolution of the normalized mass of the biggest 
MST fragment (dashed line). 

From this figure we can see that the time at which 
Ak attains its asymptotic value ( panels d, e, and f) 
signals the time at which, in average, the fragmentation 
pattern is settled for each event of a given energy and 
initial density values. To support this idea, in the upper 



panel of Figure 13 we show MST spectra calculated at 
different times. It can be noticed that a,t t ^ AOto a 
significative proportion of the mid-size asymptotic MST- 
fragments are already produced, i.e. internal surfaces 
already appeared. For later times the big residue evolves 
just by light particles evaporation. 

As a result of the expansion, physical observables like 
Ak change in time during a relatively short lapse af- 
ter which they attain almost constant values. A freeze- 
out time, Tfo, can then be introduced, after of which 
the evolution becomes almost irrelevant. Moreover, Fig- 
ures Eld, Oe, and Of also indicate that asymptotic mea- 
surements allow to reconstruct quantities as early as 
Tfo ^ 40to- 

This observation agrees with previous results. In refer- 
ence ^245 it was shown that for unconstrained systems a 
time of fragment formation (t/^), related to the stabiliza- 
tion of microscopic clusters composition, can be defined. 
For the cases studied we found that Tff ^ Tfo, meaning 
that the information is frozen once the fragment surfaces 
have developed. From the lower panels we see that the 
variance of the kinetic energy, ctk, also attains its asymp- 
totic value at this time. It is important to stress that this 
definition of a freeze out time corresponds to a MaxEnt 
principle within the constraints imposed by the dynami- 
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cal evolution. As such, this definition does not imply that 
the freeze out configuration corresponds to an absolute 
entropy maximum, i.e. to a thermodynamic equilibrium. 

For the lowest considered energy, ax does not change 
significantly in time and this stays true if we change the 
initial density. This is the expected behavior as the sys- 
tem, in almost every realization, only evaporates light 
clusters in its early evolution, and a big fragment is al- 
ways present almost independent of the system volume. 
On the other hand, for the more energetic cases the ini- 
tial decrease of ax can be related to the preequilibrium 
dynamics that involves a flux development (expansion) 
and prompt emission of light particles. Afterwards, the 
development of internal surfaces gives rise to a large set 
of possible decaying channels (opening of phase space) 
that is responsible for the uk enhancement. 

Ak asymptotic values characterize the number of dif- 
ferent fragmentation configurations compatible with the 
total energy constraint. Large kinetic energy fluctuations 
reflect large potential energy fluctuations that can be di- 
rectly related to event by event final cluster distributions. 
In this sense K fluctuations can be associated with the 
available number of decaying channels for a given energy, 
i.e. with the microcanonical entropy. 




energy (e^^) 



V. DISCUSSION 

We summarize in Fig0]the results of Ak calculations 
obtained so far. In the upper panel we show the fluc- 
tuations in kinetic energy for the three confined states 
analyzed at densities pi (circles), p2 (squares), and pa 
(diamonds). The solid lines are not a fit, but an Ak es- 
timation using information from the respective system's 
caloric curve. The nice agreement between the lines and 
the symbols demonstrates the quality of our numerical 
sampling. In the middle and lower panels we report the 
same quantity calculated at the time of uk stabilization 
(as was mentioned above, this time is well-correlated with 
the time of fragment formation) and for the correspond- 
ing asymptotic configurations respectively. In the upper 
panel the horizontal dotted-dashed line depicts the ex- 
pected canonical value for Ak- We included it also in 
the middle and bottom panels, only as a reference value. 
In the fluctuation analysis performed on experimental 
heavy ion data, the average kinetic energy is approxi- 
mately corrected for evaporation and for the Coulomb 
boost, while fluctuations are assumed to be correctly 
given by the asymptotic measurement This means 
that the middle panel is the most relevant for a compar- 
ison with experimental data. 

The confined state with density pi, that belongs to 
region A in the phase diagram of FigQ (full circles) dis- 
plays a maximum that is well above the canonical value. 
As such there is a region for which the thermal response 
function is negative. On the other hand for the confined 
states p2 and p3, that belong to region B and C respec- 
tively of the phase diagram, the corresponding values 



FIG. 4: Top: Relative kinetic energy fiuctua- 
tions as a function of the total energy for densi- 
ties pi (circles), p2 (squares), and ps (diamonds). 
Middle: same as top, but evaluated at the time 
of gk stabilization (t(pi) ~ 5io, t{P2) ^ 30io, 
^(pa) ~ 40io)- Bottom: same as top, but evalu- 
ated for asymptotic configurations. See text for details. 



of the kinetic energy fluctuations are always below the 
canonical value. 

If we now turn our attention to the results of the same 
calculation but now performed on the states resulting 
from removing the confining walls and allowing the sys- 
tem to expand and fragment, the results displayed in the 
middle and lower panel of Fig0| are obtained. Almost 
independent of the initial condition, the normalized fiuc- 
tuations fall on a curve which is mainly determined by 
the total energy and fairly close to the thermal response 
of an initially diluted system in microcanonical equilib- 
rium. The kinetic energy fluctuations that were below 
the threshold for all energies have been clearly enhanced 
in the energy region deflned by < < O.Seq. This 
feature, as was already mentioned, reflects the existence 
of quite different available decaying channels that are ex- 
plored by the systems dynamics. In particular, at these 
energies, partitions with big residues on one hand, and 
high multiplicity partitions with no such big structures 
on the other, can be produced giving rise to power-law 
like distributions. If we look at the asymptotic stage 
(lower panel), the picture is somewhat blurred by sec- 
ondary evaporation that continuously decreases the av- 
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erage kinetic energy without affecting th fiuctuation. It 
is worth noticing that at time of fragment formation the 
average system volume is approximately independent of 
the initial density: this means that for all initial condi- 
tions the information is frozen when the system reaches 
the region of the phase diagram (region A of figure Q]) 
where abnormal fluctuations are expected in equilibrium, 
i.e. in the case of a full opening of the accessible phase 
space. This observation may suggest that the collective 
motion leads the system to a low density freeze out con- 
figuration which is not far away from a thermodynamic 
equilibrium. However, to make any conclusion about the 
possible degree of equilibration and the characteristics of 
the relevant statistical ensemble, a more detailed analysis 
is needed. 



VI. CONCLUSION 

In conclusion, in this work we have shown that for 
free to expand fragmenting systems, initially equilibrated 
within sharp boundary conditions, an intermediate stage 
exists where new microscopic correlations arise. This 
happens in the presence of a collective expansive motion 
that induce internal surfaces to appear in the system, 
preceding the formation of well defined asymptotic clus- 



ters. 

In addition, the analysis of the dynamical evolution 
of Ak turned out to be fruitful. On one hand it al- 
lowed us to give a physically meaningful estimation of the 
freeze-out time, which turned out to be in good agree- 
ment with the corresponding fragment formation time. 
On the other, it was shown that a link can be established 
between asymptotic values of that observable and mass 
distribution functions. In isolated systems, like the one 
studied here, a maximum of Ak as a function of the sys- 
tem energy can be related to a noticeably enlargement of 
the number of possible decaying channels. This enhance- 
ment of the relative kinetic energy fluctuation can thus 
be associated with an opening of phase space that is ex- 
pected in phase transitions taking place in finite isolated 
systems. 

It is worth noting that even if the non-constrained frag- 
ment formation stage seems to provide a quite effective 
memory loosing mechanism (see Fig. EJ, preliminary re- 
sults support the idea that information on the initial sys- 
tem density can still be retrieved from certain asymptotic 
cluster configurational features. This analysis, along with 
a thorough characterization of the system equilibration 
properties at r/o is currently under progress. 

This work was done with partial financial supports 
from UBA and CONICET. 
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